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ABSTRACT 

We report the results of an analysis of two XMM-Newton/EPlC-pn spectra of the bright ultraluminous X-ray 
source M81 X-9 (Holmberg IX X-l), obtained in snapshot observations. Soft thermal emission is clearly revealed 
in spectra dominated by hard power-law components. Depending on the model used, M81 X-9 was observed at a 
luminosity of Lx = 1 .0- 1 .6 x 10 4l) erg s 1 (0.3-10.0keV). The variability previously observed in this source signals 
that it is an accreting source which likely harbors a black hole. Remarkably, accretion disk models for the soft 
thermal emission yield very low inner disk temperatures (kT = 0.17-0.29 keV, including 90% confidence errors 
and variations between observations and disk models), and improve the fit statistic over any single-component 
continuum model at the 6a level of confidence. This represents much stronger evidence for a cool disk than 
prior evidence which combined spectra from different observatories, and the strongest evidence of a cool disk 
in an ultraluminous X-ray source apart from NGC 1313 X-l. Like NGC 1313 X-l, scaling the temperatures 
measured in M8 1 X-9 to those commonly seen in stellar-mass Galactic black holes at their highest observed fluxes 
(kT ~ 1 keV) may imply that M81 X-9 harbors a black hole with a mass on the order of 10 3 M ; the measured 
disk component normalization and broad-band luminosity imply black hole masses on the order of 10 2 M Q . It is 
therefore possible that these sources harbor 10 3 M Q black holes accreting at Lx — 0. 1 x LecIcL.- While these results 
do not represent proof that M81 X-9 harbors an intermediate-mass black hole, radio and optical observations 
suggest that beaming and anisotropic emission from a stellar-mass black hole are unlikely to account for the 
implied luminosity. We further argue that the strength of the hard emission in these sources and well-established 
phenomena frequently observed in stellar-mass black holes near to the Eddington limit suggest that optically-thick 
photospheres are unlikely to be the origin of the cool thermal emission in bright ultraluminous X-ray sources. For 
comparison to M8 1 X-9, we have also analyzed the previously unpublished EPIC-pn spectrum of NGC 1 3 1 3 X- 1 ; 
cool disk emission is again observed and refined spectral fit parameters and mass estimates are reported. 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are point-like X-ray 
sources in nearby galaxies for which the implied luminosity ex- 
ceeds the Eddington limit (L Edd , = 1.3 x 10 38 (M /Mq) erg s" 1 ) 
for an isotropically-emitting black hole of 10 Mq (as per 
dynamically-constrained Galactic black holes; see McClintock 
& Remillard 2003). A number of these sources were first iden- 
tified with Einstein (see, e.g., Fabbiano 1989), but the spa- 
tial resolution and sensitivity of Chandra has shown that most 
ULXs are likely point sources. Short-term and longer-term 
variability studies have shown that most ULXs must be accret- 
ing sources, likely harboring black holes feeding from compan- 
ion stars in binary systems (for reviews, see Fabbiano & White 
2003; Miller & Colbert 2003). 

Intermediate-mass black holes (IMBHs, 10 2-5 Mq) provide 
an attractive explanation for ULXs. However, this interpreta- 
tion requires excellent evidence that alternative explanations 
are unlikely to hold. Viable models for ULXs which only 
require stellar-mass black holes (~ 10 Mq) and which may 
explain at least the lower-luminosity ULX sources include: 
anisotropic emission due to a funnel-like inner disk geome- 
try (King et al. 2001), strongly beamed emission due to a line 



of sight coincident with a jet axis (Reynolds 1997; Kording, 
Falcke, & Markoff 2002), and "slim" disks (Watarai, Mizuno, 
& Mineshige 2001) and/or radiation pressure-dominated super- 
Eddington accretion disks (e.g., Begelman 2002). 

While some prior analyses of X-ray spectra (e.g., with ASCA) 
found some evidence for IMBHs, the sensitivity of these ob- 
servations only statistically required single-component spec- 
tral models (see Colbert & Mushotzky 1999, Makishima et 
al. 2000), which prevented stronger conclusions. The high 
effective area and sensitivity of XMM-Newton are changing 
this observational constraint (Miller et al. 2003, Strohmayer & 
Mushotzky 2003). Similarly, high-quality optical observations 
of ULXs are now being made (e.g., Pakull & Mirioni 2002) 
which reveal relatively symmetric nebulae around some sources 
with spectra indicating the ULX is acting on the local environ- 
ment. The symmetry of these nebulae may indicate that these 
ULXs emit isotropically, and that funneling does not hold. Sim- 
ilarly, radio to X-ray flux ratios may be used to constrain rela- 
tivistic beaming models. Miller et al. (2003) analyzed XMM- 
Afewfon/EPIC-MOS spectra of the ULXs NGC 1313 X-l and 
X-2. These ULX spectra were the first to statistically require 
thermal disk and power-law continuum components. Remark- 
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ably, the disk temperatures measured in those spectra are 5- 
10 times lower than those commonly measured in stellar-mass 
black holes accreting at high rates. Those temperatures, the 
normalization of the disk components, and the very high lu- 
minosities implied suggest that NGC 1313 X-l and X-2 may 
harbor IMBHs. This suggestion is strengthened by symmet- 
ric optical nebulae around X-l and X-2 which suggest isotropic 
emission, and radio observations which likely rule-out relativis- 
tic beaming (Miller et al. 2003). XMM-Newton observations of 
the brightest ULX in M82 reveal quasi-periodic oscillations and 
even an Fe Ka emission line; the detection of these disk signa- 
tures strongly rules-out beaming and likely also funneling, and 
also represents strong evidence that the M82 ULX may harbor 
an IMBH (Strohmayer & Mushotzky 2003). 

The ULX M81 X-9 (actually in the dwarf galaxy Holm- 
berg IX, but hereafter referred to as M8 1 X-9) is an exception- 
ally bright (L x > 10 40 erg s _1 ) source which is variable on the 
timescale of weeks and months, with weak indications for a 
cool disk component (La Parola et al. 2001). The source po- 
sition is roughly 2' from the optical center of Ho IX (Paturel 
& Petit 2002). In order to better understand the nature of M81 
X-9, we analyzed spectra obtained in two short XMM-Newton 
observations available in the public archive. Herein, we re- 
port the results of fits to the XMM- Newton/EPlC-pn spectra of 
M8 1 X-9. For comparison, we also analyzed the EPIC-pn spec- 
trum of NGC 1 3 1 3 X- 1 (only the EPIC-MOS spectra were ana- 
lyzed previously). Fits to these X-ray spectra suggest that both 
sources may harbor IMBHs. We discuss these results within the 
context of other recent ULX observations, and prevalent models 
for ULX accretion flow geometries. 

2. OBSERVATIONS AND DATA REDUCTION 

In this work, we consider the time-averaged XMM- 
Newton/EPlC-pn spectra gathered in two observations of M8 1 
X-9 and a single observation of NGC 1313 X-l. The data sets 
were obtained through the XMM-Newton public data archive. 
The EPIC-pn camera has a higher effective area than the EPIC- 
MOS cameras, and drives the results of any joint spectral anal- 
ysis. The results of an analysis of the EPIC-MOS spectra of 
NGC 1313 X-l and X-2 are reported in Miller et al. (2003). 
That analysis used the XMM-Newton reduction and analysis 
suite SAS version 5.3.3, which failed to produce EPIC-pn event 
lists for the observation of NGC 1313. The reduction and anal- 
ysis reported in this work used SAS version 5.4.2, which suc- 
cessfully produced an EPIC-pn event list for the NGC 1313 ob- 
servation. The results of fits to the EPIC-pn spectrum of NGC 
1313 X-l — the stronger IMBH candidate of the two ULXs 
considered Miller et al. (2003) — are reported in this work. 

M81 X-9 (Ho IX X-l) was observed on 10 April 2002 start- 
ing at 17:37:52 (UT), and again on 16 April 2002, starting 
at 17:33:15 (UT). In both cases, the camera was operated in 
"PrimeFullWindow" mode, with the "thin" optical blocking fil- 
ter. Application of the standard time filtering resulted in net 
exposures of 7.8 ksec and 8.5 ksec, respectively. In both cases, 
events were extracted in a circle within 24" of the source posi- 
tion (09 /! 57 m 54 s ,69°03'46", La Parola et al. 2001). Due to the 
proximity of the source to the edge of the CCD and the rela- 
tively few bright sources apparent in surrounding region, back- 
ground counts were extracted in an adjacent circle with a radius 
of 24". NGC 1313 was observed by XMM-Newton on 17 Octo- 
ber 2000 starting at 03:59:23 (UT). The camera was operated in 
"PrimeFullWindow" mode with the "medium" optical blocking 



filter. Application of the standard time filtering resulted in a net 
exposure of 29.3 ksec. Events were extracted in a circle within 
24" of the position reported by Miller et al. (2003). A back- 
ground spectrum was extracted in an annulus between 24-30", 
again centered on the source position. 

Each of the spectra were made by applying the selection cri- 
teria described in the MPE "cookbook". Briefly, these selec- 
tions are as follows: we set "FLAG=0" to reject bad pixels 
and events too close to chip edges, event patterns 0^1 were al- 
lowed (these patterns correspond to "singles" and "doubles"), 
and the spectral channels were grouped by a factor of 5. The 
"canned" response files appropriate to the CCD, patterns, and 
optical blocking filter were used to fit the data. Using the SAS 
tool "epatplot" and simulating XMM-Newton/EPlC-pn spectra 
using the HEASARC tool "PIMMS" (based on prior reported 
fluxes and spectral parameters), we found that the effects of 
photon pile-up are negligible in the spectra. Prior to fitting, the 
source and background spectra were rebinned (using "grppha" 
within LHEASOFT version 5.2) to require at least 15 counts 
per bin to ensure the validity of x 2 statistical analysis. Spectral 
fits were made using XSPEC version 1 1.2 (Arnaud & Dorman 
2000). All spectral fits were made in the 0.3-10.0 keV band. 
All errors reported in this work are 90% confidence errors, and 
were evaluated while allowing all spectral model components 
to vary to take account of dependencies between parameters 
and components. For all spectral fits, the standard F test was 
employed to calculate the significance of adding model compo- 
nents. 

3. ANALYSIS AND RESULTS 

We considered a number of continuum models commonly 
applied to Galactic X-ray binary systems, modified by ab- 
sorption in the interstellar medium using the XSPEC model 
"phabs". Although somewhat technical, it is worth noting that 
fits to XMM-Newton/RGS spectra of bright Galactic X-ray bi- 
naries show that "phabs" places the neutral oxygen absorp- 
tion edge at an energy slightly below that expected for neu- 
tral atomic oxygen (0.543 keV). At CCD resolution, this can 
have the effect of creating a false emission line feature at an 
energy just above that which "phabs" assumes for the oxygen 
edge. Although different fixes for this problem are possible, 
we adopted the following approach (in part because it is easily 
reproducible). Each spectrum was fit with a model consisting 
of a simple power-law modified by "phabs"; the edge depth 
expected for the neutral oxygen edge based on this column and 
assuming solar abundances was calculated; all continuum mod- 
els then fit to the spectra were modified by "vphabs" with the 
abundance of oxygen set to zero, with an additional edge fixed 
at the calculated depth added at 0.543 keV. Regardless of the 
continuum model, this absorption model proves to be a better 
fit to each spectrum in the region of the oxygen edge, and re- 
moves the (false) apparent emission line. 

The results of all spectral fits are reported in Table 1. The 
spectra of M81 X-9 are shown in Figures 1, 2, 4, and 5, and the 
spectrum of NGC 1 3 1 3 X- 1 is shown in Figures 3 and 6. While 
a simple power-law proves to be the best single-component con- 
tinuum model in each case, it is only an acceptable model for 
the latter observation of M81 X-9. The canonical multi-color 
disk (MCD) blackbody model (Mitsuda et al. 1984) — an ap- 
proximation to the Shakura & Sunyaev (1973) disk model — is 
not an acceptable fit to any of the spectra when it is the only 
continuum component. Similarly, single-component thermal 
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bremsstrahlung models are not an acceptable fit to any of the 
spectra. Though not listed in Table 1, single-component fits 
with any diffuse thermal plasma model are not acceptable; this 
is expected given the lack of emission lines in any spectrum. 

Following Miller et al. (2003), we proceeded to consider two- 
component continuum models: MCD plus power-law, "diskpn" 
plus power-law, and "MEKAL" plus power-law. While the 
MCD model is convenient in that 1) it is a common compo- 
nent in published fits to Galactic black hole binary spectra, 2) 
its temperature and normalization allow scaling from which 
constraints on the mass of the accretor in ULXs may be esti- 
mated, it is an approximation to the Shakura & Sunyaev (1973) 
disk in that it does not include the zero-torque inner bound- 
ary condition. The "diskpn" model is a disk model that in- 
cludes corrections (based on a pseudo-Newtonian potential) for 
the temperature behavior near to the black hole (Gierliski et al. 
1999). As with the MCD model, the temperature and normal- 
ization parameters in this model permit scalings which can be 
used to constrain the mass of the accretor in ULXs. In addi- 
tion, the inner disk radius is also a free parameter in the diskpn 
model (R in = 6 R g was fixed in all fits, where R g = GM/c 2 ). 
The "MEKAL" model describes diffuse thermal plasmas; the 
results of fits with this model are not listed in Table 1 as the den- 
sity parameter of this model is totally unconstrained (likely, due 
to the total lack of credible emission lines in the spectra), and 
because "MEKAL" plus power-law models prove to be signifi- 
cantly worse than disk plus power-law models. This is particu- 
larly clear in the case of NGC 1313 X-l: the best-fit "MEKAL" 
plus power-law model gives x 2 l v = 499.0/486 (where v is the 
number of degrees of freedom); this is significantly worse than 
fits with a model consisting of an MCD plus power-law com- 
ponents (x 2 /v = 459.8/487) or diskpn plus power-law com- 
ponents (x 2 /v = 461.0/487). Finally, it is worth noting that 
archival observations suggest that NGC 1313 X-l and M81 X- 
9 are not diffuse sources, but are constistent with point sources 
at Chandra resolution (after accounting for pile-up). 

Cool disk components are formally statistically required to fit 
the first spectrum of M8 1 X-9 (the addition of a disk component 
is significant at the 6cr level of confidence) and the spectrum of 
NGC 1313 X-l (the addition of a disk component is signifi- 
cant at more than the 8cr level of confidence). Figure 2 clearly 
shows that there is also excess soft emission in the second spec- 
trum of M8 1 X-9. It was shown previously that soft excesses of 
this kind are robust against under-abundances in the local ab- 
sorbing material (Miller et al. 2003a), and that result holds for 
these spectra as well. When a disk component is included in 
the model, the parameters of the disk component can be con- 
strained and the addition gives an improvement which is signif- 
icant at nearly the la level. Remarkably, the disk components 
in fits to M81 X-9 yield very low temperatures (ranging be- 
tween kT = 0.17-0.29 keV, including 90% confidence errors, 
and considering each of the two models fit to the two observa- 
tions of M81 X-9). Inner disk temperatures ranging between 
kT = 0.21-0.25 keV are measured for the disk component in 
the spectrum of NGC 1313 X-l. The power-law indices mea- 
sured are harder than those commonly measured in the "very 
high state" in stellar-mass Galactic black holes, but softer than 
those commonly measured in the "low/hard" state. 

Estimates for the mass of the black holes in M81 X-9 and 
NGC 1313 X-l are given in Table 1, based on three indepen- 
dent methods: 1) scaling the inner disk temperature to values 
commonly measured in Galactic black holes, 2) direct estimates 



from the normalization of the disk components, and 3) scaling 
from the Eddington limit equation. 

At the highest mass accretion rates observed in stellar-mass 
black holes, inner disk temperatures near to or above kT = 
1 keV are often observed (for a review, see McClintock & 
Remillard 2003; for a recent example, see Park et al. 2003). In 
some cases, however, the inner disk temperatures can be as high 
as kT ~ 2 keV (for examples and a comparison to ULXs, see 
Makishima et al. 2000). Using the fact that in blackbody disk 
models, T oc M~ l l A , it is possible to scale temperature to mass 
via the relation (M ULX /M 1Q m ) <* (Wo M e /kT ULX ) 4 . Assum- 
ing that M8 1 X-9 and NGC 1 3 1 3 X- 1 are accreting at or near to 
the Eddington limit (the regime in which hot disks are observed 
in stellar-mass Galactic black holes), this scaling suggests that 
these ULXs harbor black holes with masses of few x 10 3 M Q 
(see Table 1). Clearly, this scaling is very sensitive to what tem- 
perature is assumed to be typical of stellar-mass Galactic black 
holes accreting at high rates. Assuming kT = 0.5 keV to be typ- 
ical gives masses of few x 10 2 M for M81 X-l and NGC 1313 
X-l; however, assuming kT = 1.0 keV to be typical for stellar- 
mass black holes accreting at or near to the Eddington limit has 
a much stronger observational basis. 

The MCD and diskpn normalizations also allow estimates of 
the black hole masses. In this case, some care is required be- 
cause it has been noted that the color temperature measured us- 
ing multi-color disk models may differ from the effective tem- 
perature (Shimura & Takahara 1995; Merloni, Fabian, & Ross 
2000). Whereas corrections divide-out in scaling the tempera- 
ture, they must be considered in scaling masses from the disk 
component normalization. The MCD normalization gives the 
inner disk extent in kilometers, which can be converted into a 
black hole mass. Following recent related work (Miller et al. 
2003; see also Sobczak et al. 2000 and Makishima et al. 2000), 
we take /?,„, corr = r\f 2 R bs. for the MCD model (where / = 1 .7 is 
the spectral hardening factor, and r\ = 0.63 is valid for i < 70° 
and accounts for the difference between the innermost radius 
and the radius of peak temperature; note also that in all disk 
normalization scalings we assumed ; = and that any higher in- 
clination gives a higher mass). For a Schwarzschild black hole, 
R in = 6 R g = 8.85 km/M Q (where R g = GM/c 1 ). Therefore, 
we can use the MCD model to obtain mass estimates via the 
relation: M ULX = r]f 2 (K/cos[i])^ 2 x (J/10 kpc) x (8.85km)" 1 
(where K is the MCD normalization, i is the inclination, and 
d is the source distance). Assuming distances of d = 3.4 Mpc 
for Ho IX (Georgiev 1991, Hill et al. 1993) and d = 3.74 Mpc 
for NGC 1313 (Tully 1998), the MCD component normaliza- 
tion suggests that M81 X-9 and NGC 1313 X-l may harbor 
black holes with masses of few x 10 2 M (see Table 1). The 
normalization of the diskpn component explicitly includes the 
mass and the color correction factor /. In estimating masses 
with the diskpn model, / = 1 .7 was again assumed, as per the 
MCD model. The diskpn normalizations measured for M81 
X-9 and NGC 1313 X-l again suggest that these ULXs may 
harbor black holes with masses of few x 100 M Q (see Table 1). 

The mass of the putative black holes powering M81 X- 
9 and NGC 1313 X-l can also be estimated by scaling 
the implied source luminosities to the Eddington limit for 
isotropically emitting accretion-powered sources (L E dd. = 1 - 3 x 
10 38 (M/M Q ) erg s _1 / M ©; Frank > Kin S> & Raine 2002). To 
make this scaling, the unabsorbed 0.3-10.0 keV flux was ex- 
trapolated to the 0.05-100.0 keV band, as the latter is more 
representative of a bolometric luminosity. This scaling again 
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implies that M81 X-9 and NGC 1313 X-l may harbor black 
holes with masses of few x 10 2 M Q (see Table 1). It is inter- 
esting to note that scaling the ULX inner disk temperatures to 
kT = 1.0 keV implies masses on the order of 10 3 M Q while 
the disk normalizations and overall luminosities imply masses 
on the order of 10 2 M Q . This may imply that these ULXs are 
10 3 M Q black holes seen at L x ~ 0. 1 x L Edd ,. 

Finally, we note that while the results obtained here for 
NGC 1313 X-l broadly agree with those previously reported 
in Miller et al. (2003), they do not agree exactly. The differ- 
ences can be attributed to several factors, including: the more 
conservative fitting range used in this analysis (0.3-10.0 keV) 
relative to the prior analysis, the higher effective area, sensitiv- 
ity, and energy resolution of the EPIC-pn camera (used in this 
analysis) relative to the EPIC-MOS cameras (used in the prior 
analysis), and the improved interstellar absorption model used 
in this analysis. The most important difference is that the disk 
temperatures measured in this analysis are 0.08 keV higher than 
in the previous analysis (kT = 0.23 keV versus kT = 0.15 keV), 
resulting in new mass estimates which are five times lower than 
prior estimates using inner disk temperatures. 

4. DISCUSSION 

We have investigated the nature of the ULX M8 1 X-9 (Holm- 
berg IX X-l) by analyzing two XMM-Newton/EPlC-pn spectra 
of the source. The most important result of this work is that cool 
thermal continuum emission is unambiguously required for the 
first time. Optically-thick disk components yield improvements 
in the fit statistic which are significant at more than the 6cr level 
of confidence. This represents the clearest evidence for a cool 
accretion disk in a ULX apart from NGC 1313 X-l (Miller et 
al. 2003; see Table 1). 

Remarkably, the inner disk temperatures measured (kT = 
0.17-0.29 keV, including 90% confidence errors on two disk 
models fit to the two observations) are well below those com- 
monly measured in stellar-mass black holes (kT ~ 1 keV, see 
McClintock & Remillard for a review). Scaling the tempera- 
tures measured implies that M8 1 X-9 may harbor a black hole 
with a mass on the order of 10 3 M Q , and scaling the disk com- 
ponent normalizations implies masses on the order of 10 2 Mq. 
It is unclear which scaling method is superior, but both scalings 
suggest that M8 1 X-9 may harbor an IMBH. Both scaling meth- 
ods can be correct if the source harbors a black hole with a mass 
on the order of 10 3 M Q which is observed at L x — 0.1 x L E dd.- 
For comparison, we also report results from an analysis of the 
EPIC-pn spectrum of NGC 1313 X- 1 , which may also harbor 
an IMBH (Miller et al. 2003). The spectra of these sources bear 
striking similarities (see Table 1). 

The long-term and short-term flux variability of M81 X-9 
demonstrated clearly by La Parola et al. (2001) establishes that 
M8 1 X-9 is an accreting source. While it is natural to ascribe 
cool thermal emission to an optically-thick accretion disk in 
accreting systems, given the implications of these cool disks it 
is especially important to understand why disk emission is the 
most viable explanation. 

Relativistic beaming — which might prevent detection of 
disk signatures — is not likely at work in M81 X-9. First, ra- 
dio observations of Ho IX by Bash & Kaufman (1986) find a 
flux density of at most 1 mJy at 20 cm at the position of the 
ULX, corresponding to a luminosity of ~ 2 x 10 34 erg s. If the 
source was active at its present level during that time (as sug- 
gested by the long-term behavior reported by La Parola et al. 



2001), this radio luminosity gives a radio to X-ray luminosity 
ratio of ~ 2 x 10~ 6 . This is at odds with a relativistic beam- 
ing scenario, as beaming should create flat spectra (see, e.g., 
Fossati et al. 1998). Second, the peak radio to X-ray flux ra- 
tios observed in Galactic X-ray binary systems are always be- 
low 2.3 x 10" 5 (Fender & Kuulkers 2001; Barth, Ho, & Sargent 
2003). The fact that the radio to X-ray flux ratio in M81 X- 
9 is likely an order of magnitude below the highest ratio seen 
in Galactic systems suggests that M8 1 X-9 is not likely to be 
a stellar-mass black hole with relativistically-beamed emission 
coincident with our line of sight to the source. 

Slightly anisotropic X-ray emission — perhaps due to a 
funnel-like geometry in the inner disk — might create a very 
hot inner disk, especially if the source of hard X-ray emission 
is central (e.g., through Comptonization). Clearly, the spectra 
we have analyzed rule-out hot disk signatures. Pakull & Miri- 
oni (2002) have reported optical nebulae around a few ULXs, 
including NGC 1313 X-l and M81 X-9. The line ratios ob- 
served in these nebulae suggest that X-ray photoionization may 
be important, and hint that some ULXs may act on their lo- 
cal ISM. The symmetry of the nebula around NGC 1313 X-l 
suggests that the source emits isotropically. While the nebula 
imaged around M8 1 X-9 is better described by an ellipse than a 
circle, the ratio of the axes is likely less than 2:1, and certainly 
inconsistent with a ratio of 10:1 ("funneled" scenarios require 
an anisotropy parameter of approximately 10; see King et al. 
2001). Thus, present data suggests that anisotropic emission is 
not the best explanation for the high flux observed from M8 1 
X-9. 

Radiation pressure-dominated accretion disks — which 
might be able to produce super-Eddington fluxes through small- 
scale photon bubble instabilities (Begelman 2002) — might be 
expected to have especially high temperatures (perhaps similar 
to the "slim" disk solution described by Watarai, Mizuno, & 
Mineshige 2001). As the spectra of M81 X-9 and NGC 1313 
X-l are inconsistent with hot disk components, it is not likely 
that such models describe these sources. 

Finally, it has recently been suggested that soft thermal 
emission in the brightest ULXs may be due to an optically- 
thick, outflowing photosphere originating at 100 Rs c hw. around 
a stellar-mass black hole accreting at a super-Eddington rate 
(e.g., King 2003). This model is inconsistent with the relative 
importance of hard X-ray emission in sources like M8 1 X-9 and 
NGC 1313 X-l. It is also inconsistent with a number of long- 
established properties observed in stellar-mass Galactic black 
holes accreting near to, at, or slightly in excess of their isotropic 
Eddington luminosities. We will briefly consider these issues 
here. 

The mechanical power in a photosphere outflowing at veloc- 
ity v should scale with the luminosity radiated by the photo- 
sphere as L mec h. ~ L ra d x (v/c), and that the photospheric radius 
is given by R p hot — (c/v) 2 x Rs c hw. (King 2003). Hard X-ray 
emission in accreting sources is usually tied to regions deep 
within the gravitational potential, as this is a convenient en- 
ergy reservoir. However, the innermost accretion regime would 
be blocked from view by an optically-thick photosphere. The 
observed hard component in such a scenario must be gener- 
ated in shocks above the photosphere, which means that the 
strength of the hard component must be driven by the outflow: 
Lhard — Lmech- Measurements indicate that the hard components 
in M81 X-9 and NGC 1313 X-l are 3-4 times stronger than 
the soft (putatively photospheric) components in each spectrum 
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(see Table 1). Even if we only assume that Lh ar d — Aso/i — L m d, 
the only way L mech ~ L rad x (v/c) holds is if (v/c) ~ 1, which 
implies that R pho , ~ (c/v) 2 x /?5 cAw . ~ 1 

Moreover, even at the highest implied accretion rates, a num- 
ber of spectral and variability properties observed in stellar- 
mass Galactic black holes would be screened if an optically- 
thick photosphere developed at R p hot — 100 Rschw.- First, hot 
(kT ~ 1 keV) thermal spectral components are observed in ev- 
ery source at peak X-ray flux; this emission can only be as- 
sociated with the inner disk as a radially distant photosphere 
should be much cooler. Second, QPOs - whether high fre- 
quency (few x 100 Hz) or low frequency (few x 1 Hz) - are 
almost certainly tied to the disk, and likely represent Keple- 
rian orbital frequencies and/or coordinate resonance frequen- 
cies. QPOs are seen over a factor of 10 3 in flux in stellar- 
mass Galactic black holes but are preferentially seen at high 
fluxes. Most importantly, QPOs are a higher fraction of the rms 
noise at high energies - QPOs are intimately tied to hard X-ray 
emission. It is very unlikely that such periodicities originate in 
shocks above an optically-thick photosphere. Third, broad, rel- 
ativistic Fe Ka emission lines are also observed in stellar-mass 
Galactic black holes over a factor of 10 3 in flux, but preferen- 
tially at the very highest fluxes (e.g., in the very high state). 
These lines are certainly tied to hard X-ray emission (see, e.g., 
Zdziarski, Lubinski, & Smith 1999), and have been observed 
to vary at frequencies as high as 6 Hz (Gilfanov, Churazov, & 
Revnivtsev 2000). These facts - and their smeared shape (very 
likely due to the extreme Doppler shifts and gravitational red- 
shift at the inner disk) - tie broad Fe Ka emission lines to the 
innermost accetion flow. The reader is directed to McClintock 
& Remillard (2003), and references therein, for a discussion of 
the stellar-mass black hole properties noted here. 

It is worth noting that partial covering models for NGC 1313 
X-l and M81 X-9 cannot be ruled-out statistically in the lim- 
ited energy range considered. Models consisting of a strong ab- 
sorber (Nh — 10 23 cm" 2 ) covering 50-70% of a steep power-law 
source (T ~ 2.7) can fit the data as well as our disk plus power- 
law models. There are several reasons, however, why partial 
covering is likely an unphysical model for these sources. First, 
partial covering is only required in Galactic black hole sources 
during X-ray "dips" (e.g., 4U 1630-472 and GRO J1655-40, 
Kuulkers et al. 1998). When partial covering is actually re- 
quired to fit the spectra of Galactic black holes and AGN, fits 
with simple models require very strong neutral Fe K absorp- 
tion edges: Ueda et al. (1998) find an edge with t = 3.3 in 
GRO J 1655-40, and Boiler et al. (2002) report an edge with 
t = 1.8 in 1H 0707-495. In contrast, X-ray dips have never 
been observed in NGC 1313 X-l or M81 X-9, and the 95% 
confidence upper limits on the strength of a neutral Fe K edges 
in the spectra of NGC 1313 X-l and M81 X-9 are r < 0.4 and 
r < 0.2, respectively. Next, it has been shown that while par- 
tial covering models can provide acceptable fits to AGN data 
below 10 keV, simultaneous high energy spectra strongly re- 
ject such models (see, e.g., Reynolds et al. 2004). Finally, the 
above discussion of photospheric models shows that an alterna- 
tive geometry wherein partial covering might arise naturally is 



very unlikely to describe sources like NGC 1 3 1 3 X- 1 and M8 1 
X-9. 

Thus, at present, the most viable explanation for the soft ther- 
mal excess observed in sources like M81 X-9 and NGC 1313 
X-l is emission from a cool accretion disk around an IMBH. 
Recent observations of bright ULXs with XMM-Newton and 
Chandra have revealed cool thermal (likely disk) components 
in a number of ULXs, including (but not limited to) NGC 1313 
X-l and X-2 (Miller et al. 2003), NGC 5408 X-l (Kaaret et al. 
2003), NGC 4038/4039 X-37 (Miller et al. 2003b), and M74 
X-l (Krauss et al. 2003). The brightest ULX in M82 may or 
may not have a cool disk; the low energy portion of the spec- 
trum from that ULX has proved difficult to measure due to pho- 
ton pile-up in CCD spectrometers (with Chandra) and diffuse 
emission in the crowded field (with XMM-Newton). However, 
beaming is confidently ruled-out in this source through the ob- 
servation of QPOs and an Fe Ka emission line (Strohmayer & 
Mushotzky 2003), and this source is an excellent IMBH can- 
didate. At most, these observations suggest that sources at 
the upper-end of the ULX luminosity distribution may harbor 
IMBHs. Many ULXs at the lower end of the luminosity dis- 
tribution may be stellar-mass black holes. As observations of 
ULXs continue to be made (especially, multi-wavelength ob- 
servations), it is likely that sub-classes of beamed, funneled, 
and IMBH sources will be distinguished. 

If M81 X-9 harbors an IMBH, we can speculate about how 
such a black hole was made. This speculation is particu- 
larly interesting because Ho IX is a dwarf galaxy. Madau & 
Rees (2001) have proposed IMBHs might have been created by 
the death of extremely massive, low-metallicity Population III 
stars. While most IMBHs created in this way may have been 
dragged to galactic centers by dynamical friction, a fraction 
might be visible today as ULXs if they have captured a donor 
star. If M81 X-9 has a mass on the order of 10 2 M Q and this 
scenario can hold in the case of dwarf galaxies, M8 1 X-9 could 
be a Population III remnant. If M8 1 X-9 has a mass on the 
order of 10 3 M , it is more likely that it grew through merg- 
ers, either in a young cluster (e.g., Ebisuzaki et al. 2001) or in 
a globular cluster (Miller & Hamilton 2002). Certainly, Miller 
(1995) notes that the vicinity of M81 X-9 may contain a num- 
ber of massive, young stars, and may show evidence for a recent 
supernova history. 

J. M. M. acknowledges support from the NSF through its As- 
tronomy and Astrophysics Postdoctoral Fellowship program, 
and useful discussions with A. Kong and P. Slane. M. C. M. 
was supported in part by NSF grant AST 00-98436, and by 
NASA grant NAG 5-13229. This work is based on observa- 
tions obtained with XMM-Newton, an ES A mission with instru- 
ments and contributions directly funded by ES A member states 
and the US (NASA). This work has made use of the tools and 
services available through HEASARC online service, which is 
operated by GSFC for NASA. The authors note for the curi- 
ous that J. M. Miller (Miller et al. 2003a, 2003b), M. C. Miller 
(Miller & Hamilton 2002, Miller & Colbert 2003), and B. F. 
Miller (Miller 1995) are indeed distinct, unrelated people. 



REFERENCES 



Arnaud, K. A., and Dorman, B., 2000, XSPEC is available via the HEASARC 

on-line service, provided by NASA/GSFC 
Barth, A. J., Ho, L. C, & Sargent, W. L. W., 2003, to appear in "AGN Physics 

with the Sloan Digital Sky Survey", eds. G. T. Richards and P. B. Hall (ASP 

Conference Series), astro-ph/03 10435 



Bash, E N., & Kaufman, M., 1986, ApJ, 310, 624 

Begelman, M., 2002, ApJ, 568, L97 

Boiler, TH., et al. 2002, MNRAS, 329, LI 

Colbert, E. J. M., & Mushotzky, R. E, 1999, ApJ, 519, 89 

Ebisuzaki, T., et al., 2001, ApJ, 562, L19 



6 



M81 X-9 (Holmberg IX X-l) 



Fabbiano, G., 1989, ARA& A, 28, 87 

Fabbiano, G., & White, N. E., 2003, to appear in "Compact Stellar X-ray 
Sources," eds. W. H. G. Lewin and M. van der Klis, Cambridge: Cambridge 
University press, astro-ph/0307142 

Fender, R. P., & Kuulkers, E., 2001, MNRAS, 324, 923 

Fossati, G., Maraschi, L., Celotti, A., Comastri, A., & Ghisellini, G., 1998, 
MNRAS, 299, 433 

Frank, J., King, A. R., & Raine, D., 2002, in 'Accretion Power in Astro- 
physics", Cambridge: Cambridge Univ. Press 

Georgiev, Ts. B., Bilkina, B. I., Tikhonov, N. A., & Karachentsev, I., 1991, A& 
AS, 89, 529 

Gierlinski, M., Zdziarski, A. A., Poutanen, J., Coppi, P. S., Ebisaa, K., & John- 
son, W. N., 1999, MNRAS, 309, 497 
Gilfanov, M., Churazov, E., & Revnivtsev, M., 2000, MNRAS, 316, 923 
Hill, J. K., et al. 1993, ApJ, 402, L45 

Kaaret, P., Corbel, S., Prestwich, A. H., & Zezas, A., 2003, Science, 299, 365 
King, A. R., 2003, Proc of the II BeppoSAX Meeting: "The Restless High En- 
ergy Universe" (Amsterdam, May 5-8 2003), eds. E. P. J. van den Heuvel, J. 
J. J. M. in 't Zand, and R. A. M. J. Wijers, astro-ph/0309524 
King, A. R., Davies, M. B., Ward, M. J., Fabbiano, G., & Elvis, M., 2001, ApJ, 
552, LI 09 

Kording, E., Falcke, H, & Markoff, S., 2002, A& A, 382, L13 
Krauss, M., et al., 2003, ApJ, submitted 

Kuulkers, E., Wijnands, R., Belloni, T., Mendez, M., van der Klis, M., & van 

Paradijs, J., 1998, ApJ, 494, 753 
La Parola, V., Peres, G., Fabbiano, G., Kim, D. W, & Bocchino, f., 2001, ApJ, 

556, 47 

Madau, P., & Rees, M. J., 2001, ApJ, 551, L27 
Makishima, K, et al. 2000, MNRAS, 313, 193 

McClintock, J. E., & Remillard, R. A., 2003, to appear in "Compact Stellar 
X-ray Sources," eds. W. H. G. Lewin and M. van der Klis, Cambridge: 
Cambridge University press, astro-ph/0306213 



Merloni, A. A., Fabian, A. C, & Ross, R. R., 2000, MNRAS, 313, 193 
Miller, B. F, 1995, ApJ, 446, 175 

Miller, J. M., Fabbiano, G., Miller, M. C, & Fabian, A. C, 2003a, ApJ, 535, 
L37 

Miller, J. M., Zezas, A., Fabbiano, G., & Schweizer, F, 2003b, ApJ, subm., 
astro-ph/0302535 

Miller, M. C, & Colbert, E. J. M., 2003, "Intermediate-Mass Black Holes", re- 
view subm to International Journal of Modern Physics D., astro-ph/0308402 
Miller, M. C, & Hamilton, D. P., 2002, MNRAS, 330, 232 
Mitsuda, K, et al., 1984, PASJ, 36, 741 

Pakull, M. W., & Mirioni, L., 2003, in "New Visions of the X-ray Universe in 
the XMM-Newton and Chandra Era" (ESA SP-488; Noordwijk: ESA), in 
press, astro-ph/0202488 

Park, S. Q., et al. 2003, ApJ, subm., astro-ph/0308363 

Paturel, G., & Petit, C, 2002, LEDA 

Reynolds, C. S., Loan, A. J., Fabian, A. C, Makishima, K, Brandt, W. N., & 

Mizuno, T., 1997, MNRAS, 286, 349 
Reynolds, C. S., Wilms, J., Begelman, M. C, Staubert, R., & Kendziorra, E., 

2004, MNRAs, subm., astro-ph/0401305 
Shakura, N. I., & Sunyaev, R. A., 1973, A& A, 86, 121 
Shimura, T., & Takahara, E, 1995, ApJ, 445, 780 
Sobczak, G. J., et al. 2000, ApJ, 544, 993 
Strohmayer, T. E., & Mushotzky, R. F, 2003, ApJ, 586, L61 
Tully, R. B., 1988, Nearby Galaxies Catalog (Cambridge: Cambridge Univ. 

Press) 

Ueda, Y., Inoue, H, Tanaka, Y., Ebisawa, K, Nagase, F, Kotani, T., & Gehrels, 

N., 1998, ApJ, 492, 782 
Watarai, K, Mizuno, T., & Mineshige, S., 2001, ApJ, 549, L77 
Zdziarski, A. A., Lubinski, P., & Smith, D. A., 1999, MNRAS, 303, Lll 



Miller, Fabian, & Miller 



7 



Table 1 

Spectral Fit Parameters 



Model/Parameter 




M81 X-9 Obs. 1 


M81 X-9 Obs. 2 


NGC 1313 X-l 


power-law 










N H (10 21 cm- 2 ) 




2.1(1) 


2.5(2) 


2.4(2) 


r 




1.93(4) 


2.00(4) 


2.04(5) 


Norm. (10~ 3 ) 




1.24(4) 


1.57(5) 


0.69(3) 


X 2 /dof 




460.7/443 


571.7/579 


555.7/489 


MCD 










N H (10 21 cm- 2 ) 




2.1 


2.5 


2.4 


kT (keV) 




0.91(5) 


1.63(4) 


0.78(5) 


Norm. (KT 1 ) 




3.5(2) 


0.40(2) 


3.2(2) 


X 2 /dof 




1823/444 


1421/579 


1919/490 


bremsstrahlung 










N H (10 21 cm- 2 ) 




2.1 


2.5 


2.4 


kT (keV) 




31(1) 


43(1) 


4.6(1) 


Norm. (1(T 3 ) 




1.04(2) 


1.22(2) 


0.70(2) 


X 2 /dof 




1510/443 


2435/579 


773/489 


MCD + power-law 










Nh (W cm z ) 




2.3(3) 


2.9(3) 


3.1(3) 


kT (keV) 




?6 


0.21(4) 


0.23(2) 


Norm. 




zu -10 


60+ 70 


28(5) 


r 




1.73(8) 


1.86(6) 


1.76(7) 


Norm. (10~ 3 ) 




0.96(9) 


1 <i ^+0 OQ 
1,J -0.05 


0.49(4) 


X 2 /dof 




419.0/441 


524.3/577 


459.8/487 


disk significance 




6.1(7 


6. So- 


> 8(7 


F (10~ 12 erg cnf 2 r 1 


) 


8+ 2 


lO! 2 


4.3(4) 


F power— law / ^tolal 




0.85 


0.84 


0.74 


io.3-lo(10 40 erg.v-') 




1 14-0 1 


1 3+0- 3 
-0.2 


0.6(1) 


£0.05-100 (10 40 erg s~ 


') 


-0-1 

2000! 3 ™ 

140™ 
330+140 
-100 
^/M-1+50 
200+,[j 


2 q+0.6 

-0.3 

„_L*7nnn 

SOOO +/UUU 
JUUU -2000 


4000 
<+ u -1000 

,,,111 no 

220+00 


M kT =i.o (Mq) 




M kT =0.5 (%) 




320+ 4 ! 


M Norm , (M©) 




J ' -240 

220$ 


400+ 4 ° 
110™ 


diskpn + power-law 










Afo (10 21 cm' 2 ) 




2.3(3) 


2.8(3) 


3.0(3) 


kT (keV) 




0.24(5) 


0.21(4) 


0.22(3) 


Norm. (1(T 4 ) 




4+16 


7+ 30 


5+ 05 


r 




i~ 74 +0.07 
1 -'^-0.04 


1.85(6) 


1.76(7) 


Norm. (10~ 3 ) 




1.0(1) 


1 3+ 02 

1 "'-0.1 


0.49(5) 


X 2 Mof 




418.8/441 


524.9/577 


461.0/487 


disk significance 




5.9(7 


6.7ct 


> 8(7 


F (10 -12 erg cm~ 2 s -1 


) 




10+ 2 


4.3(4) 


F power-law jFtotal 




0.85 


0.85 


0.74 


L .3-io(10 40 erg.r 1 ) 




1 1+ 03 


1 3+ ' 3 

'-0.2 


0.6(1) 


£0.05-100 (10 40 erg s~ 


l ) 


9 7+0.7 

z -0.3 


7 q+0.6 
-0.3, 


1 4+0.2 
400o|°°° 


A%- = i. (M ) 




3000^™ 


5000+7000 
-2000 


A*w=o.5 (M Q ) 




200 -i«o 
330+ 28 
-20 


300+ 4 * 
-100 


300+ 2 <J 
240+10 


Mjv orm . (M Q ) 




380™ 
2201 


m l x Il em . (Mq) 




200™ 


1 10+^ 



NOTE. — The results of fitting simple models to the EPIC-pn spectra of M81 X-9 and NGC 1313 X-l. The XSPEC model "vphabs"was used to measure the 
equivalent neutral hydrogen column density along the line of sight. All errors reported above are 90% confidence errors. The significance of adding disk components 
was calculated using the standard F test. Errors in parentheses are the error in the last digit. Fluxes are "unabsorbed" fluxes, and the reported luminosities assume 
a distance of 3.4 Mpc for Ho IX and 3.7 Mpc for NGC 1313. Mass estimates based on scaling the inner disk temperatures to those commonly measured in 10 Mq 

-1/4 

black holes accreting at high rates assume T oc M gH . Scalings are given assuming typical inner disk temperatures of kT = 1.0 keV and kT = 0.5 keV, although 
the literature favors the former. Please see the text for details on mass estimates obtained through the disk component normalizations. Mass estimates based on the 
Eddington luminosity make use of the 0.05-100 keV luminosity (extrapolated from the 0.3-10.0 keV flux) and assume L Edd = 1.3 X 10 38 (M/Mq) erg/s). 
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M81 X-9 (Holmberg IX X-l) 



in 
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FlG. 1 .— The ratio of the EPIC-pn spectrum of M81 X-9 obtained in the first observation to a model consisting of only power-law continuum fit in the 3.0-10.0 keV 
band. The soft continuum excess is readily apparent at low energy. The spectrum shown has been rebinned for visual clarity. 

CM 
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FlG. 2. — The ratio of the EPIC-pn spectrum of M81 X-9 obtained in the second observation to a model consisting of only power-law continuum fit in the 
3.0-10.0 keV band. The soft continuum excess is readily apparent at low energy. The spectrum shown has been rebinned for visual clarity. 
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FlG. 3. — The ratio of the EPIC-pn spectrum of NGC 1313 X-l obtained in the second observation to a model consisting of only power-law continuum fit in the 
3.0-10.0 keV band. The soft continuum excess is readily apparent at low energy. The spectrum shown has been rebinned for visual clarity. 
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FlG. 4. — The EPIC-pn spectrum of M81 X-9 obtained with the first observation of the source. The spectrum has been fit with a model consisting of MCD (blue) 
and power-law (red) continuum components, and the data/model ratio is shown below. The shape of the disk component and ratio plot are nearly identical if the 
"diskpn" disk model is fit instead of the MCD model. See Table 1 for details of the spectral model. 



|_! j \ 

1 * 

-r 

i i 


^X 

. 1 — J — ■ 











1 2 
Energy (keV) 



FlG. 5. — The EPIC-pn spectrum of M81 X-9 obtained with the second observation of the source. The spectrum has been fit with a model consisting of MCD 
(blue) and power-law (red) continuum components, and the data/model ratio is shown below. The shape of the disk component and ratio plot are nearly identical if 
the "diskpn" disk model is lit instead of the MCD model. See Table 1 for details of the spectral model. 




FlG. 6. — The EPIC-pn spectrum of NGC 1313 X-l obtained with the first observation of the source. The spectrum has been fit with a model consisting of MCD 
(blue) and power-law (red) continuum components, and the data/model ratio is shown below. The shape of the disk component and ratio plot are nearly identical if 
the "diskpn" disk model is fit instead of the MCD model. See Table 1 for details of the spectral model. 



